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Abstract 


We study thermally induced density perturbations during inflation. This 
scenario is characterized by two thermodynamical conditions: (1) The pri¬ 
mordial perturbations originate in the epoch when the inflationary universe 
contains a thermalized heat bath. (2) The perturbations of the inflationary 
scalar field are given by the fluctuation-dissipation relation. We show that 
the spectrum of the primordial perturbations is of power law, but tilted, and 
there is a relation between the amplitude and the index of the power spec¬ 
trum. Aside from the mass scale of the inflation, the amplitude-index relation 
does not depend on other parameters like (^-factor. These results are found to 
be well consistent with observations of the temperature fluctuations of cosmic 
microwave background if the mass scale of the inflation is about 10^® GeV. 
Instead of the purely adiabatic case, the consequent density perturbation is 
an admixture of adiabatic and isocurvature one. Therefore, the detection 
of super-Hubble suppression of the spectrum would be effective for further 

discrimination between the thermally originated models and others. 
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I. INTRODUCTION 


Recently, the possibility of thermally induced initial density perturbations in an inflation¬ 
ary cosmology was examined |^. It was found that for a wide parameter range within the 
standard inflation model, the thermal fluctuations of the scalar held can play an important 
and even dominant role in producing the initial perturbations. 

The “standard” scenario of inhation assumes that during the inhation the cosmic 
expansion was dominated by the vacuum energy of a scalar held 0, and the thermalized 
component has never existed, or if a thermal component like radiation existed at the onset 
of inhation, it was blown oh by the exponential expansion. The thermal component was not 
established or re-established until heating or reheating caused by the decay and thermaliza- 
tion of the 0 held into relativistic particles. The thermal inhation scenario also assumes that 
the inhation is dominated by a scalar held (j). Unlike the standard model, the thermal sce¬ 
nario requires that a thermalized component be established or re-established before the end 
of the inhation. Therefore, the thermal model can be considered similar to the “standard” 
model, but with a heating or reheating stage of the 0-decayed particles while the inhation 
was still in progress. 

It has been realized that if the interaction of the inhation held with light particles is weak, 
namely, of the order of gravitational strength or slightly stronger, the thermal component 
cannot maintain itself during the inhation. In contrast, if the interaction is strong enough 
to lead to thermalization, it is possible to maintain a heat bath in the epoch of inhation 
P]. Therefore, the thermal scenario probably can be “ installed” onto various versions of 
inhation by adding a stronger interaction between the 0 held and the heat bath. 

In this paper, we study the properties of density huctuations of thermal origin. These 
density perturbations are characterized by two thermodynamical conditions: (1) the pertur¬ 
bations originate in the epoch when the inhationary universe contains a thermalized heat 
bath; (2) the perturbations of the inhationary scalar held are given by the huctuation- 
dissipation relation. Therefore, one can expect that, like many problems of thermal physics. 
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some features of the perturbations should be mainly determined by thermodynamical re¬ 
quirements, but less dependent on the details of the model, such as the inflationary potential. 
Obviously, these features would be useful in testing the thermal scenario of perturbation pro¬ 
duction, and discriminating it from other models. 

Up to now, the only possible way of directly testing models of the origin of cosmic 
inhomogeneity is to compare the predicted spectrum of the perturbations with the observed 
temperature fluctuations of the cosmic microwave background (CMB). Therefore, we will 
focus on the observable features of the perturbations, including the power law index, the 
magnitude of the spectrum, and the behavior of the perturbations on scales larger than the 
Hubble radius. 

In Sec. II we discuss the exact and approximate solutions of the evolution of the radiation 
component in the thermal inflation scenario. Sec. Ill gives the calculation of the features of 
the density perturbations produced during the thermal inflation. Finally, the observational 
tests and possible discrimination between this model and others are discussed in Sec. IV. 

II. HEATING OR REHEATING DURING INFLATION 

A. Thermal inflation scenario 

The dynamics of a universe consisting of an inflationary scalar held 0 and a thermal- 
ized radiation component can be described by the following three equations |]^. First, the 
equation of motion of the 0 held is 

0-f3i/0 + F^0 + U'(0) = 0 . (2.1) 

where V (0) is the effective potential. Second, the equation of the radiation component (heat 
bath) is given by the hrst law of thermodynamics as 

Pr + 4:H pr = T , (2.2) 

where pr is the energy density of the thermal component. The temperature of the thermal 
bath can be calculated by pr = (7r^/30)5fefr7"‘^, Qes being the effective number of degrees of 
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freedom at temperature T. The friction terms of T,^ in Eqs. (2.1) and (2.2) phenomenolog¬ 
ically describe the interaction between 0 held and the heat bath (radiation component) [Q. 
The third equation is the Hubble “constant” H = R/R, i.e, 

+ ly + r(^)]. (2.3) 

where mpi = 1/a/G is the Planck mass. 

In the standard scenario of inhation, the radiation component is assumed to be nonex¬ 
istent, or to be blown away by the exponential expansion, and the universe is cool (T H) 
during the inhation era. To establish or re-establish a radiation component, a heating or 
reheating epoch following the inhation is needed. However, if the friction coefficient T^ is 
comparable to or larger than H, and there is strong coupling between the decay products of 
the scalar held, then the system of the decaying products will quickly reach thermal equi¬ 
librium. In this case, a non-zero and thermal component pr can be maintained during the 
inhation. 

One may be concerned that the friction term is not appropriate for describing the 
energy transfer via 0 held particle production under far-from-equilibrium conditions. Indeed, 
it has been shown that in some models of reheating, the hrst stage of the 0 held decaying 
into other bosons is due to a parametric resonance |^. However, in the last stage while the 
thermalization is established the friction mechanism prevails. Therefore, the description of 
T^ is reasonable if there is a thermalized heat bath |Q]. 

Briehy, thermal inhation scenario can be characterized by the following properties: 

1. ) The relaxing time scale of the heat bath is shorter than the expansion 

> H, (2.4) 

2. ) The energy of the universe is dominated by 0 held potential 

1/(0) » 0V2, H(0) » Pr. (2.5) 

3. ) The temperature of the heat bath is larger than H, 
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T> H. 


( 2 . 6 ) 


4.) As in the standard inflation scenario, Eq. (2.1) should have a slow-roll solution. That 
is, the potential of the inflaton cj) held is “flat”, 

|l/'(0)|«(3if + r^)E'/2^0), (2.7) 

and 


|V'"(.^)|«(3ff + r^)2. 


(2.8) 


As such, the slow-roll solution of Eq. (2.1) is given by 

i., v'Cf) 


(2,9) 


3H -|- 

It has been shown that conditions (2.4) - (2.8) can hold simultaneously if the coupling 
constant T^ and the potential V (0) of inflation field satisfy |l| 


( M M / M \M 

Vmpi/ W Vmpi/ IE’ 


( 2 . 10 ) 


< E'(0) < mplV{(j)). (2.11) 

where M is the mass scale of the inflation, and IE = 0^/2E(0) is the ratio of the kinetic 
and potential energy of the (j) held. In the case of (M/mpi)^ <C 1, there is a large region in 
the parameter space of the potential fulfilling the conditions Eqs. (2.10) and (2.11). 

Exactly speaking, we should use a finite temperature effective potential V{(j), T) in the 
dynamics of thermal inflation. However, for many popular potentials, the correction due to 
finite temperature is negligible. As an example, we consider a (j)'^ potential given by 

E(0) = A(02 - ^2)2. (2.12) 


The slow-roll solution gives 



Stt V (0) 
3 rripi 


(2.13) 
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where V^(0) = = Xa^, and the subscript i denotes the time of the beginning of the 

inflation. The condition (2.8) can certainly be satished if A < (M/mpi)^. On the other 
hand, the leading temperature correction of the potential (2.12) is AT^^^. Therefore, AT^ < 

H^, i.e, the influence of the hnite temperature effective potential 

is negligible when 0 < mpp 

B. Evolution of Radiation Component 

The evolution of the held in the thermal inhation scenario is about the same as in 
the standard model. However, the radiation component shows very different behavior in 
these two scenarios. In the standard model, the radiation component did not exist before 
reheating, while in the thermal scenario, it appeared quite early. To illuminate this point, we 
hrst plot a typical numerical solution of Eqs. (2.1) - (2.3) in Fig. 1, where the potential V (0) 
is given by (2.12) and the relevant parameters are taken to be M = 10^^ GeV, A = 4 ■ 10“^^, 

= ISJIi and ges = 100. It should be pointed out that g^eff-factor generally is a function of 
T, and the unknown function ge{{(T) will lead to uncertainty in the solutions. Fortunately, 
the variation of g^eff-factor has only a slight effect on the problems we are going to study. 

The solution shown in Fig. 1 is easy to understand. As expected, the inflation is 
dominated by the vacuum energy of the (j) held, and H remains almost constant during the 
inhation. The evolution of the radiation contains three phases. Phase 1 covers the period 
during which the radiation temperature T drops drastically due to the inhationary expansion. 
This phase is the same as in cool inhation models. When t = tb, T will rebound, and the 
evolution of radiation evolves into phase 2 which lasts from tb to tf. In this phase, H 
while the radiation temperature increases due to the friction of the 0 held. Therefore, the 
second phase is actually the epoch of inhation plus reheating. Finally, phase 3 begins at 
t = tf when the energy of the radiation components is comparable to the potential energy of 
0 held, and the inhationary process stops. The third phase is, in fact, the standard radiation 
dominated universe, and radiation temperature follows the standard evolution of adiabatic 
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decreases with the expansion of the universe. 

Therefore, the only difference between the thermal and the standard scenarios is that 
the “reheating” starts at the time tb < tf when the inflation still prevails. In other words, 
the reheating stage is merged into the epoch of inflation. A specihc reheating stage between 
the epochs of inflation and radiation dominated universe is not needed. The transition from 
inflation to a radiation dominated universe is now smooth, because a heat bath has already 
been established before the end of the inflationary expansion. 

Moreover, the solution of phase 2 is independent of the initial conditions of the radiation 
and the 0 held. Whatever the initial condition of the radiation, the solution during phase 2 
is always about the same, as shown in Fig. 1, because the radiation in phase 2 is produced 
during inhation, and the initial radiation has almost been consumed by the expansion. 
Therefore, this solution of phase 2 can be installed in both initially hot and cool models. It 
can be used as either a reheating or heating solution. 

The features shown in Fig. 1 are typical of thermal inhation. One can describe these 
features by approximate solutions of Eqs. (2.1)-(2.3). During the slow-roll stage we have 
0 (T, it would be reasonable to neglect the 0^ term in V'{(j)). Eq. (2.1) is then 

0 + (3hf + F^)0 - 4Aa20 = 0. (2.14) 

The inhationary solution of Eq. (2.14) can immediately be found as 

0 = 0*e“^*, (2.15) 

where 0* is the initial value of 0 held. 0* can be estimated by E'(00 = Hf, which gives 
(pi ~ (M/mpi)^(T. The uncertainty of pi slightly ahects the results discussed below. 

The coefficient a in Eq. (2.15) is given by 

where 7 = (F,^ + 3H)/Hi. Since we have required A < {M/m-p\Y (Af/'^rpi)^, and 7 3> 1, 

Eq. (2.16) gives a ~ 3A^/^(mpi/M)^/27r7 1. 
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From Eq. (2.15), one can find the exact solution of Eq. (2.2) as 


Pr 




(2.17) 


where C = a‘^HTip(j)f/2{a + 2). If the initial condition is taken to be Pr(0) = aM^, where 
a = (7r^/30)5feff , we have C = Pr(0) — C. It is clear from Eq. (2.17) that when the hrst term 
on the right-hand side of Eq. (2.17) is dominant, the evolution is in phase 1; and phase 2 
starts when the second term becomes dominant. The transfer from phase 1 to phase 2 occurs 
at the rebound of the temperature of radiation. The transition time ti, can be determined 
by {dp/dt)tf, = 0. We have 


Htb ~ 


1 

2(a + 2) 


In 


/4(a + 2) aM^\ 
a^H 


(2.18) 


The temperature at the rebound is 


Tb = a 




(2.19) 


Therefore, the radiation temperature T in phases 1 (t < tb) and 2 (t > tb) can be 
approximately expressed as 


T(t) 


TbC if t <tb 

jf t > tb 


( 2 . 20 ) 


Substituting Eqs. (2.15), (2.18), and (2.19) into the approximate solution (2.20) of t > U, 
we have 


= aT^ = 

^ 4 ^ AH ^ AH^ 


( 2 . 21 ) 


Substituting (2.21) into Eq. (2.2), we have 


Pr ~ 0. 


( 2 . 22 ) 


Solution (2.22) implies that the evolution of the radiation is simply determined by the 
balance between the depletion of the initial radiation density due to inflationary expansion 
and the production of thermal components by the friction. 



From Eq. (2.21), we have T = The second phase will end at 

the time tf, when 0 approaches the minimum of V{(j)), i.e, 0 < (t/\/2. Accordingly, the 
temperature at the end of phase 2 should be 

Tf < = eM, (2.23) 


where 

1/4 

( 2 . 

Eq. (2.23) shows that the temperature of radiation will rise to of order M when the inflation 
ends. Accordingly, both the heating (or reheating) and the inflation are underway in phase 
2, and no post-inflation reheating is needed. 

Using the same parameters as the numerical solution of Fig. 1, the approximate solutions 
(2.20) and (2.23) are plotted. It shows that this approximation is in excellent agreement 
with the numerical one. 



ciF(^ 

2a(F0 -1- 3H) 


C. Duration of the thermal inflation 

Generally, the rebound temperature can be either greater or less than H. Figure 2 
plots a typical solution in the case of Tf, < H. This solution shows that the temperature T 
is lower than H in the period of tf, < t < te, and it is higher than H, i.e, T > H when t > te- 
From Fig. 2, the evolution of T at Tf, < shows about the same behavior as that in the 
period of Tf, > H. However, we should not consider the solution to be physical when T < H, 
because it is impossible to dehne a thermalized heat bath with temperature less than the 
Hawking temperature H. Nevertheless, this solution should be available when t > and 
T> H. 

As has been discussed in Sec. H.B, the evolution of radiation at t > te is independent of 
whether the radiation existed at t < te- The behavior of T at the period t > te is completely 
determined by the competition between the diluting and producing radiation at t > tf,. 
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Initial information about the radiation has been washed out during the thermalization. For 
instance, one can replace the solution at t < tg by a standard (cool) inflation. Therefore, 
te can be considered as the time of the onset of thermal inflation, and the duration of the 
thermal inflation is — 4- The number of e-folds of the growth of the cosmic factor R 
during the thermal inflation is then 


N = 


b 2 Tt 

Hdt ~ — In —. 
a H 


(2.25) 


One can also formally calculate the number of e-folds of the growth in phase 2 as 

ftf 2 Tf 

N2= Hdt ~ -In;^, 

Jtb O' l}j 


(2.26) 


and the number of e-folds of the total growth as 


Nt = 


, 2, /eM\ 2 

Hdt ~ - In + Hth = — 


a 


6 = - In 
a 


1/4 


M 


-) ^ 


(2.27) 


As expected, both N 2 and Nt depend on the initial value of the held but N does not. 
Therefore, 0* will not lead to uncertainty in our analysis if we only study the problems of 
thermal evolution at t^, < t < tf. 


III. THE PRIMORDIAL DENSITY PERTURBATIONS 

A. Fluctuations of the (j) field 

The possibility for the fiuctuations of thermal origin during infiation is simply because 
there is a gap between the energy scales of the infiation M, and the Hawking temperature 
H. There is room for the temperature of the radiation component T to be less than M 
but larger than H. The condition T < M ensures the condition Eq. (2.5), and then, the 
evolution of the 0 field leads to an infiation. Condition T > if is necessary in order to have 
the perturbations be dominated by thermal fiuctuations. 

To calculate the thermal fiuctuations, we rewrite Eq. (2.9) as 

d4 _ 1 dF[0] 

dt 3H -1- T^ d(j) ’ 
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(3.1) 







where F[0] = V^(0). This is, in fact, a rate equation of order parameter 0, describing a 
homogeneous system with the free energy F[0] approaching to equilibrium. A rate equation 
like (3.1) cannot correctly describe the approach to equilibrium during a phase transition 
without a noise term |^. To ensure that the system approaches the global minimum, the 
order parameter dynamics is not pure relaxation but exhibits fluctuations, arising from the 
microscopic degrees of freedom. Namely, Eq. (3.1) should have an additional noise term r] 
such as 


dfj) 1 SF 


(3.2) 


This is the Langevin equation for a system with one degree of freedom, and the correlation 
functions of r] is determined by the fluctuation-dissipation theorem as 


{v{t)v{t')) 


3H^ T 
271 3H -|- r<^ 


(3.3) 


The fluctuations of 0 held for E"(0) <C are then 


m)^} = ^HT. 


(3.4) 


Obviously, when T > H, the thermal huctuations are larger than the quantum ones. When 
T < H, the thermal component is negligible, and huctuations should be of quantum origin, 

z.e, mr) = 

Like all inhation models, the primordial perturbations undergo an evolution of crossing 
outside, and later again crossing back inside the Hubble radius. The density perturbations 
at re-entering the horizon are 

(L) = , (3.5) 

\p)h 0^ + (4/3)p^ 

All quantities in the right-hand side of Eq. (3.5) are calculated at the time when the relevant 
perturbations crossed out beyond the horizon at the inhationary epoch. 

We are only interested in the primordial density perturbations produced in the period 
of thermal inhation, i.e, t > t^. From Eq. (3.5), we have 
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(3.6) 


\p),^ 2\2it^) \g,s-P) ' 


where the temperature T should be in the range Tf > T > H. 

B. Power law index 

Since the inflation is immediately followed by the radiation dominated epoch, with no 
reheating in between, the comoving scale of a perturbation with crossing over (the Hubble 
radius) at time t is given by 


— = 27r—— 

Ho Ho Tf 


(3.7) 


where Tq and Ho are the present CMB temperature and Hubble constant, respectively. Eq. 
(3.7) shows that the smaller t is, the smaller k will be. This is the so-called “first out - last 
in” of the evolution of density perturbations produced by the inflation [§]. 

Using solution (2.20) and (3.7), the perturbations (3.6) can be rewritten as 




(3.8) 


where k^ is the wavenumber of perturbations crossing out of horizon at tg. It is 


Tf 


Tf 


(3.9) 


Therefore, the primordial density perturbations produced during thermal inflation are of 
power law with an index —3a/2. While in terms of the density perturbations at a given 
time t, it can also be expressed as a power law as 0 

/(^pV 


oc 


P 


k^+^, ii k> ke, 


(3.10) 


and then the index n becomes 


n = 1 — 3 q(/2. 


(3.11) 


Consequently, the power spectrum of thermally induced perturbations is a tilted power law. 
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The thermal scenario requires that all perturbations re-entering in the Hubble radius 
originate in the period of thermal inflation. The longest wavelength of the perturbation 
(3.8), i.e, 2ti/ kei should be larger than the present Hubble radius. In other words, N should 
be large enough in order to have the inflated patch cover the entire Hubble radius. From 
Eq. (3.9), we have then 


N >\n 


' HTq ' 

HoTf, 


In^ ~ 70. 
-no 


(3.12) 


where we used (Tq/Hq) 3> (Tf/H). Thus, from Eqs. (2.25) and (3.11), we can hud a lower 
limit to n as 


^min 1 3 


HTf/H) 


(3.13) 


HTo/HoY 

Figure 3 plots rirmn as a function of the inflation mass scale M and Urmn is sensitive to 
M, but only slightly varying when ^feff is in the range 10^ — 2 ■ 10^. Thus, as mentioned in 
Sec. H.B, the uncertainty given by g^efr-factor is not important. 

We now try to hnd an upper limit to n. From Eqs. (2.9) and (2.25), and ~ 
87rH(0)/3m|i, we have 

N= j*' f*‘ Yd4>, (3.14) 

(f) 3mpi V 

where 0e = 0* exp(Q;i7te) and 0/ = 0* exp(Q;i7t/) ~ For the potential (2.12), Eq. 

(3.14) gives 


N = 


271'y , a 


Srripi JiPf 


= 2 


Try 

37TTpj 


4 2 


Because 0e ^ Eq- (3.15) can be rewritten as 


(3.15) 


N=l^j{4aN-l) 

b mi 


6pi 


(3.16) 


or 


a = 


1 -|- 6iVm|i/7r(T^7 

AN 


> 


/^pi\ 


2 3 

> 


27r7 V a / 27r7’ 


(3.17) 
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here we used < m|[. Hence, an upper limit to n is found to be 

n < 1 - (3.18) 

47r7 

For most interesting cases, 7 > 70 (see Sec. III.C), and therefore, n < 0.99. Eq. (3.17) also 
provides another lower limit, a > l/4iV, and the corresponding upper limit to n is 

n<l-3/8iV. (3.19) 

In the case of thermal inflation, iV > 70 [Eq. (3.12)], and therefore,we obtain the same 
upper limit as Eq. (3.18), i.e, n < 0.99. Thus, roughly speaking, the allowed values of the 
power law index should fall within the area between the curves of rimm and the line n = 1 
in Fig. 3. 


C. Amplitudes of perturbations 


To calculate the amplitude of the perturbations we rewrite spectrum (3.8) into 


'5p 


k 


n—1 


p ) ]h 


, ii k > ke, 


(3.20) 


where ko = A is the normalized amplitnde of the pertnrbation on scale k = ko, which 

re-enters the Hubble radius 1/Hq at present time. From Eqs. (3.6), and (3.7), we have 


A = 


405 F^ / HqTj 


H 




(3.21) 


(27rr g,„H \ HTo J \T. 

Applying Eq. (2.20), the radiation temperatnre at the horizon-crossing moment, t, can be 
expressed as T{t) = T/exp[(n — — t)/3]. With the help of Eq. (2.25), we obtain 


^ 0 ^ 

(27r)3^eff77 VToA \h) 

_ 405 F 0 / 


A = 


(27r)3^eff77 VTo. 

On the other hand, Eqs. (2.23), (2.24), and (2.25) render 


(3.22) 


9a(f 


H (i-n)-3a(f)4e4(i-W3' 


(3.23) 
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Substituting Eq. (3.23) into Eq. (3.22), we have finally 


A = 


405 Ott^ /Ho 


n—1 


' H \A^(l—n)nN/3 




(3.24) 


(27r)3 30 V To ; (1 - n) - 3 a(f )%4(i-n)iv/3' 

Because for a given energy scale of inflation M, power law index n and Qes, N can be 
approximately determined by Eq. (2.25), the amplitude A is a function of M, n and 
Figure 4 shows A as a function of n for given parameters M and Qes- The values of n in Fig. 
4 cover the entire allowed range shown in Fig. 3. Again, it can be seen from Fig. 4 that 
gives only a slight effect on the magnitude of A. But the magnitude of A will change by 
about 8 orders when n changes from 1 to 0.60. 

To check the consistency of these solutions, we calculate T^ as function of M, n and 
The result is plotted in Fig. 5. Figure 5 shows that in all cases T^ will be in the range 
around lO^hf. These values are satisfied very well the condition of FT < T^ < M. Figure 5 
also shows that in the parameter range allowed by the thermal scenario the values of T,^ are 
mainly determined by M, but not n and ges- 

To replace n, A and T^ can also be treated as functions of M, N, and ges- Figures 6 and 
7 give the relationships of A - and T,^ - A^, respectively. Similarly, the variation of ges 
does not have much effect on the amplitudes of primordial density perturbations, and the 
friction T,*. 


D. Suppression on snper-Hubble scales 

As in the standard scenario, the fluctuations of (j) field in the thermal scenario leads 
to curvature perturbations. However, the radiation in the thermal scenario will also 
contribute to the total perturbations. Instead of the adiabatic mode, the fluctuation in pr 
will provide an isocurvatnre component of pertnrbations. This is because the relaxation 
time-scale of the radiation during the thermal inflation is eqnal to or less than l/T,^, and 
mnch shorter than 1/H. Namely, the radiation component is always in thermal eqnilibrium, 
and their thermal (temperatnre) fiuctnation is much less than the fiuctnations caused by 0 
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field phase transition (flnctnation-dissipation relation). Therefore, in terms of the eqnation 
of state, the componnd system of (j) held and radiation is non-nniform. This is isocnrvatnre 
initial condition which specihes the relative abnndance of the (/)-decayed particles and the 
thermal component. 

One can simply estimate the amplitnde of the isocnrvatnre pertnrbations by the following 


way |TI[. As has been discnssed in Sec. II, the inhation ended at the stage of pr ~ p</,. Hence, 
at f/ when the radiation dominated phase starts, the total density contrast shonld be 


^ ^ + 5 Pr ^ 1 5p^ 

P P(f) T Pr 2 p^p 

The amplitnde of the entropy pertnrbations is then given by 


(3.25) 


- 


3 5p</, 


3 ?,5p^ 3 6p 

-dr — - - — - 


(3.26) 


4 4 S Pr 4 p 

These isocnrvatnre hnctnations will eventnally lead to density pertnrbations of amplitndes 
S once they enter the horizon. 

Accordingly, the amplitnde of density pertnrbations before re-entering the horizon shonld 
be lower than the re-entered pertnrbations by a factor of abont 1 -|- 3/4 = 1.75. In other 
words, the density pertnrbations are snppressed on snper-Hnbble scales. The power spec- 
trnm shonld therefore be modihed to 


( h \ 


(3.27) 


where the fnnction f{k) is 


fik) = 


where s = 1 — (1.75) 


-1 


1 k > ko 

1 — s k < ko 

0.4 is called the snper-Hnbble snppression factor. 


(3.28) 


IV. DISCUSSIONS AND CONCLUSIONS 


Assnming that the initial inhomogeneity of the nniverse originated in the period of ther¬ 
mal inflation, the properties of the mass density pertnrbations have been stndied. A nnmber 
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of the features of the perturbations have been accounted for, including (1) the spectrum of 
the perturbations is of power law, but titled, and the allowed values of the power law in¬ 
dex are in the range of about 0.6 < n < 1; (2) for a given mass scale of the inflation, the 
amplitude of the perturbations A is almost completely determined by the index n; (3) the 
perturbations are suppressed on super-Hubble scales. 

In this analysis the popular 0^ potential of the inflationary held was employed. However, 
only one parameter, the mass scale of the inhation M, is involved in the predictions listed 
above, and no other details of the model are needed. Therefore, we believe that these 
qualitative features should be common for the thermal inhation scenario. Actually, the 
thermal model is based on two thermodynamical requirements: (a) a thermalized heat bath 
exists during inhation and (b) the initial huctuations are given by the huctuation-dissipation 
theorem. Like many thermodynamical problems, the results are mainly determined by the 
thermodynamical conditions, and less dependent on other details. 

The predictions (1) and (2) can be tested by current observations of the temperature 
huctuations of the Comics Microwave Background (CMB). The four-year COBE-DMR ob¬ 
servations found n ~ 1.2 ± 0.3, Qrms-ps ~ or Qrms-ps ~ ISpiC for hxed n = 1 

H- We have plotted this observational allowed area of (A, n) in Fig. 4. Meanwhile, the 
relation between n and Qrms-PS, he. 5Q/Q = —(ln5)5n from data htting, has also been 
plotted in Fig. 4. All these results are consistent with the thermal scenario if the mass scale 
of the inhation is equal to about 10^^ GeV. 

However, one should be careful in the comparison of the thermal inhation predictions 
with the COBE results. In all the htting of the COBE-DMR observations done so far, only 
the spectra described by two parameters (n, Qrms-Ps) or (n. A) has been considered. In 
other words, the published results of the COBE data analysis provide only the most likely 
area in the bi-parameter (n, Qrms-Ps) space. This htting implies that it was assumed that 
there is no suppression on scales larger than the Hubble radius, or s = 0 in Eqs. (3.27) and 
(3.28). 

To test the predictions of thermal inhation, we need to ht the data with spectrum 
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suppressed on the super-Hubble scales. The most likely values of n and Qrms-ps (or n and 
A) are generally dependent on the suppression on super-Hubble scales. Because long wave 
contributions are reduced by the super-Hubble suppression, one can expect that the most 
likely value of n given by htting with the suppression will be smaller than that without the 
suppression []K|. The result would then be moving toward the direction in favor of a tilted 
spectrum. Nevertheless, the suppression should not dramatically change the allowed area 
shown in Fig. 4. Therefore, one may already conclude that for thermal inflation, the most 
likely value of the inflationary mass scale is about 10^® GeV, and the regions larger than 
10^® GeV or lower than 10^^ GeV are likely to be excluded. 

For a better test, we need a tri-parameter £t to hnd the most likely values of n. A, and s. 
Various models of the perturbation origins are difficult to distinguish in the bi-parameter n 
and Qrms-PS (^) space. All models of the “standard” inflation, thermal inflation, and phase 
transition defects, can produce the perturbation spectrum with about same bi-parameters. 
However, the super-Hubble suppression factor s is different for different models. In the 
“standard” inflation model, the perturbations are assumed to be adiabatic, and then s = 0. 
On the other hand, isocurvature perturbations require s = 1. For models such as defects in 
phase transitions, causality allows no perturbations on scales larger than the Hubble radius, 
we have s = 1 too. The thermally originated fluctuations are neither purely adiabatic, nor 
purely isocurvature. Because the thermal inflation ends at the time when the energy of 
the radiation is comparable to the 0 held energy, the suppression factor s should not be 
restricted to 0 or 1. Therefore, the value of s is crucial to distinguish the thermal model 
from the others. 
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FIGURES 

FIG. 1. Numerical solutions of the evolutions of (j) field and radiations in a thermal inflation 
model with U(<^) = \{4>^ — at the inflation mass scale M = 10^® GeV. The parameters are 
taken to be A = 4 • 10“^^, T,^ = ISRj, and = 100. The dot-dashed and solid lines are the 
evolutions of H{t) and the temperature T of the radiation, respectively. The circles are given 
by the approximate solution of T, which perfectly matches the exact solution by the end of the 
slow-roll era. tb is the time at which the temperature rebounds, and tf is the end of inflation. T 
and H are in units of Hi = [87rU(0)/3mpj]^/^, and t is in units of 1/Hi. 

FIG. 2. An exact solution of H and T similar to that shown in Figure 1, but the parameters are 
taken to be A = 4 • 10“^'^, = 150Rj, = 100. te is the time of the onset of thermal inflation. 

FIG. 3. The lower limit of the spectral index n as a function of the mass scale M. g^g is taken 
to be 100 (dotted line), 500 (dot-dashed line) and 2000 (solid line), respectively. 

FIG. 4. The amplitudes of the power spectrum as a function of n which covers the allowed 
range nmin < n < 1. The mass scales M are taken to be 10^^, 10^^ and 10^® GeV labeled at the 
three sets of curves, and are 100 (dotted lines), 500 (dot-dashed lines) and 2000 (solid lines). 
The shaded zone is the allowed area of the values (n and A) given by the GOBE-DMR data. The 
solid line within the shaded zone represents the approximated relation 5Q/Q = — (ln5)(In. 

FIG. 5. T,^ as a function of spectral index n, which covers the allowed range nmin < n < 1. 
The mass scales M are taken to be 10^^, 10^^ and 10^® GeV labeled at the three sets of curves, 
and ges are 100 (dotted lines), 500 (dot-dashed lines) and 2000 (solid lines). T^ is in units of Hi. 

FIG. 6. The amplitude of the power spectrum as a function of the thermal duration N in the 
range of V > ln(To/Ro) ~ 70, i.e, the range of thermal inflation. The mass scales M are taken to 
be 10^^, 10^^ and 10^® GeV labeled at the three sets of curves, and g^Q are 100 (dotted lines), 500 
(dot-dashed lines) and 2000 (solid lines). 
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FIG. 7. as a function of the thermal duration N in the range of > ln(To/i7o) ~ 70. The 
mass scales M are taken to be 10^^, 10^^ and 10^® GeV labeled at the three sets of curves, and ges 
are 100 (dotted lines), 500 (dot-dashed lines) and 2000 (solid lines). F^ is in units of Hi. 
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